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Abstract 

The concept of a very long baseline neutrino experiment with quasi monochro- 
matic neutrino beam and very large area underground detector is discussed. The 
detector could be placed in the existing 20 km tunnel at IHEP, Protvino. The 
High Intensity Proton Accelerators (HIPA) which are planned to be built in Japan 
(JAEPJ-KEK, baseline of ~ 7000 km) and Germany (GSI, baseline of ~ 2000 km) as 
well as the Main Injector at Fermilab (~ 7600 km) are considered as possible sources 
of neutrino beams. The oscillations are analysed in the three-neutrino scheme tak- 
ing into account terrestrial matter effects. In the proposed experiment it is feasible 
to observe the oscillation pattern as an unique proof of the existence of neutrino 
oscillations. Precise measurements of disappearance oscillation parameters of the 
muon neutrinos and antineutrinos can be done within a reasonable time. 



1 Introduction 



The observation of the atmospheric disappearance fl], [| stimulated an impressive num- 
ber of proposals of new experiments with neutrino beams from proton synchrotrons [|3], [|, 
HI H 0! H an d from neutrino factories f|. Main goals of this new generation of experi- 
ments are the ultimate confirmation of neutrino oscillations and precise measurements of 
oscillation parameters. 

Very long baseline (VLBL) neutrino-oscillation experiments with very large area (VLA) 
detectors are especially effective to observe the neutrino oscillation patterns and to mea- 
sure precisely oscillation parameters. Here we consider the prospects to use underground 
detector placed in the existing UNK tunnel to carry out these measurements. As a source 
for neutrino beams we consider the high intensity proton accelerators which are planned to 
be built in Japan (JAERI-KEK, Tokaimura, baseline of ~ 7000 km) and Germany (GSI, 
Darmstadt, baseline of ~ 2000 km) as well as the existing Main Injector at Fermilab 
(~ 7600 km). 

In the second section the physics justification of the proposed approach is given. The 
concept of a possible experimental lay-out including a neutrino focusing system and UNK 
underground detector are described in the third section. Physics performance is outlined 
in the fourth section. The last section is devoted to conclusion. 

2 Physics motivation 

We propose to measure the and survival probabilities P(z/ M (z/ M ) — > v^ip^)) as a 
function of the neutrino energy. The UNK underground detector will register muons 
from and charged current (CC) interactions in surrounding soil mainly. Expected 
effective mass will be at a ~ 1 MTon level for a few GeV neutrino energy. 

The motivation of the oscillation approach is based on the observation of atmospheric 
Vfj, disappearance and a small level, if any, of reactor v e disappearance. Existing exper- 
imental data for the muon (electron) neutrino disappearance are analyzed in terms of 
the simple formula for oscillation probability P(z/ At (^ e ) — > v^{v e )) = 1 — 1^ ■ sin 2 (1.27 • 
Am 2 ■ L/E) between two neutrino of different flavours (see Table 1 for the experimental 
restrictions on the parameters). 



Table 1. Summary of experimental results on neutrino disappearance. 



Experiments 


Best values 


Allowed region 


disappearance of 
atmospheric |10 


J„=l. 
Am 2 = 2.5 • 10" 3 eV 2 


J M > 0.84 (99% CL) 
1.2 • 10- 3 < Am 2 < 4 • 10- 3 eV 2 (99% CL) 


disappearance of 
reactor u e [[□]] 




I e < 0.1 (90% CL) 



At Am 2 = 2.5 • 10~ 3 eV 2 and with a few GeV neutrino energy the oscillation length 
is in the region of 10 3 — 10 4 km. Therefore, the measurement of the P(i/ m (i/ m ) — > /^(^v)) 
probability as a function of neutrino energy is a natural way to observe the neutrino 
oscillation pattern and to measure precisely the oscillation parameters for this pattern in 
the VLBL neutrino experiments. 
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The observation of atmospheric disappearance is just an indication on neutrino 
oscillations but not a definite proof of it. Within the proposed approach different aspects 
of neutrino physics could be studied: 

• Measurement of and disappearance patterns. 

In the case of a simple oscillation pattern Piy^ — > v^) = l — I^-sm 2 (1.27-Am 2 -L/E) 
the precise measurements of Am 2 and 1^ can be done for and v^. 

• Search for no oscillations or non-standard oscillations. 

The models of neutrino disappearance by decay [|1^] or nonstandard oscillations due 
to flavour changing neutrino interactions could be checked with high sensitivity. 
Observation of more complicated oscillation pattern could be a signal, for example, 
of large extra dimensions [|14[] . 

• Measurement of the difference A M = P(z/ M — > v^) — Pip^ — > z/ M ). 

It can come from a fake CPT violation due to the terrestrial matter effects [15| and 
can be from a genuine CPT violation if it exists |[L6|| . As fake CPT violation depends 
on the CP-violating phase, the proposed experiment could provide information on 
it as well. 

Below we consider the simple oscillation model in vacuum as the baseline, concentrat- 
ing on the ability to observe the oscillation pattern. 

Fig. 1 shows the and v e survival probabilities as functions of L/E. Arrows indicate 
the maximal L/E values which can be reached with the MINOS |J, the CERN-Gran 
Sasso (CNGS) [|, the JHF-SK and the CERN Neutrino Factory (CNF) [| planned 
experiments. 
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Fig. 1. Piv^ — > Vp) and P[y e — > v e ) survival probabilities. 

As it is seen in the figure CNF @, CNGS § and MINOS § projects are at the 
beginning of the first wave of the oscillation curve, where oscillation patterns are not 
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seen. The JHF-SK project |7j] has some chance to observe an oscillation pattern for 
Am 2 > 3 • 10~ 3 eV 2 . 

As one can see in the figure the oscillation pattern in principle could be measured for 
the atmospheric v^. However to do this, a precise reconstruction of the neutrino direction 
is needed to obtain the neutrino pass length. In the SK and MACRO detectors the neu- 
trino direction is defined from the direction of muons which have wide angular distribution 
at low neutrino energy ({E„ ) = 2.4 GeV). Therefore, only the average disappearance was 
observed. Thus the oscillation curve is clearly visible if 

• oscillation phase <f> osc = 1.27 • Am 2 • L/E ~ (2n + 1) • tt/2, n = 0, 1, 

• error S(p osc < 7r/4. 

Translation of these conditions to the energy resolution gives 

5E V < 1 

E v < 2(2n + l)' 

It means that the soft restriction is for n = and 1, i.e., for the first and the second 
minima. 

If we look at the Earth globes, the 1-st condition can be fulfilled if neutrino beams are 
sent to Protvino from the Fermilab MI, the JHF or the GSI (see expected JHF-UNK and 
GSI-UNK L/E regions on Fig. 1). In the Table 2 the possible characteristics of beams 
from the MI, the JHF and the GSI are presented. 

Table 2. Comparison of neutrino sources for UNK VLBL neutrino oscillation experiment. 



Parameter 


JHF, Japan 


GSI, Germany* 


MI FNAL, USA 


Baseline from the UNK, km 


7000 


2000 


7600 


Tilt angle, degree 


33 


9 


34 


Proton momentum, GeV/c 


50 


50 


120 


Cycle time, s 


3 


3 


2 


One turn time, [is 


5 


4 


10 


Proton intensity/spill 


3.3 ■ 10 i4 


1 ■ 10 i4 


4 ■ 10 i3 


Protons/s 


1 ■ 10 14 


3 • 10 13 


2 ■ 10 13 


7r + yield/p at 7 GeV/c 


0.05 


0.05 


0.12 


7T+ yield/s at 7 GeV/c 


5 • 10 12 


1.7- 10 12 


2.4 • 10 12 


duty factor 


1.7- 10~ 6 


1.3 • 10" 6 


5 • 10~ 6 



* assumed parameters 



Below we consider the JHF and the GSI as neutrino beam sites in comparison. 
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3 Experimental lay-out 



3.1 Concept of experiment 

Sketch of the general experimental lay-out is shown in Fig. 2. 
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Fig. 2. General experimental lay-out. 

Far neutrino source from the JHF (baseline of 7000 km) or GSI (baseline of 2000 km) 
have to point to the UNK tunnel which is ~ 50 m underground in average. For such kind 
of distances, to have enough statistics, we should think about a detector with effective 
mass at the 1 MTon level. In such a case it is possible to use the surrounding soil of 
the UNK tunnel as neutrino target. It is proposed to cover the UNK tunnel walls by 
scintillation counters. For such an experimental concept, the UNK detector, which can 
be considered as a huge double scintillator counter, will count mainly muons from iv^) 
CC interactions in the surrounding soil. 

The circumference length of the UNK tunnel is ~ 21 km, and the diameter is ~ 5 m. 
It means that an area of about ~ 10 5 m 2 can be used for a very long baseline neutrino 
experiment. Dimensions and time resolution of individual scintillation counters will allow 
us to select muons from the expected neutrino direction. 

To suppress the cosmic ray background the fast (one turn) extraction of the proton 
beam and the synchronization of a neutrino source and of the UNK detector [|17|] are 
ultimately needed. 

To observe the neutrino oscillation pattern we intend to use a few (at least three) 
neutrino energy settings with a reasonable energy resolution and just count the double 
coincidences in the UNK detector. Ratios of the measured to expected counts at different 
energy settings will provide the possibility to observe and oscillation patterns in the 
disappearance mode. 

If we assume that at least the n = 1 minimum should be detected then 

SB " < 0.15. 



E. 



V 



The ideal disappearance curves and also the smeared curves for SE U /E U = 0.36 (a WB 
like beam) and 5E U /E U = 0.15 (a NB like beam) as a function of the neutrino energy for 
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Am 2 = 2.5 ■ 10 3 eV 2 for the GSI and JHF cases are shown in Fig. 3. 




Fig. 3. Ideal (solid line), smeared for SE U /E U = 0.36 (dashed line) and for 5E U /E U = 0.15 
(dotted line) oscillation curves for the GSI and the JHF cases. 

The importance of the SE U /E U = 0.15 condition for a clear observation of the neutrino 
oscillation pattern is seen from this figure. 

The proposed concept can be considered as a further development of the atmospheric 
neutrino experiments with upward- down-going muons. Here we already know: 

• the exact distance from source to detector; 

• the type of neutrino [y^ or v^)\ 

• the neutrino energy. 



3.2 Concept of the neutrino focusing system 

To estimate the possible and event rates in the UNK detector we used the focusing 
system based on the magnetic horns, which were optimized for the NuMI Project 
This focusing system (Fig. 4) can provide desired SE U /E l/ = 0.15 condition. 

The 50 GeV primary proton beam is used to produce pions in the ~ 2 interaction 
lengths graphite target T. Two 3 m long magnetic horns H1-H2 with parabolic shaped 
inner conductors are used to focus the resulting pion beam down a drift space where pions 
will decay to muon neutrinos. Four 2 m long dipoles B1-B2 with 400 mm gaps are used 
to obtain a beam of secondaries with relatively small Ap/p. The total length of the decay 
region is equal to 400 m including the 40 m length target area and the 360 m length decay 
pipe of 1 m radius. 

This focusing system may be tuned to different neutrino energy ranges by scaling 
the dipole currents and by corresponding adjustment of the target location, keeping the 
current in both horns at its nominal value of 200 kA. Energy spectra of and CC 
events in the UNK detector are shown in Fig. 5 for a few different neutrino energy settings. 
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Fig. 4. The layout of the focusing system. 




Fig. 5. Energy spectra of the (solid lines) and (dashed lines) CC events in the 
UNK detector for the 1.7, 3.5, 7.0, 10.5 and 14.0 GeV energy settings. 



Table 3. Focusing system and neutrino beam parameters. Target positions are given 
with respect to the upstream end of the first horn. The baseline is 7000 km. 



Beam energy range, GeV 


1.0-2.3 


2.3-5.0 


5.0-9.0 


6.0-15. 


8.0-19. 


Target Z M p S ^ ream , m 


-0.34 


-0.34 


-1.94 


-2.64 


-3.34 


Mean energy, GeV 


1.7 


3.5 


7.0 


10.5 


14.0 


Energy spread HWHH, % 


~16 


~16 


~14 


~14 


~13 


events/kTon/10 21 pot 
Fraction of P M events, % 


0.45 


1.5 


2.8 


2.1 


0.33 


0.8 


0.4 


0.3 


<0.3 


<0.3 


events/kTon/10 21 pot 
Fraction of events, % 


0.11 


0.38 


0.49 


0.27 


0.033 


<16 


<6 


<10 


<17 


<30 
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The optics layout, event rates and backgrounds are summarized in Table 3. 
One should note the following: 

1. Three tune-ups of the considered focusing system span together the energy range 
from 1.7 to 7 GeV and tune-ups intermediate between these three can be easily 
achieved. The more complicated configuration of a focusing system (with obvi- 
ously smaller focusing efficiency for some particular tune-up) is required to provide 
neutrino beams with (E u ) from 1.7 up to 13-14 GeV. 

2. To obtain more intensive neutrino beams with ~ 3 and 6 events/kTon/10 21 pot for 
(E u ) = 3.5 and 7.0 GeV, respectively, one may use a rectilinear focusing system with 
horns at the same positions but without dipoles. In this case the energy spreads of 
beams (HWHH) are of about 37%. 

3.3 Concept of detector 

The proposed baseline design for the UNK detector relies on the known plastic scintillator 
extrusion technology. The UNK detector uses extruded plastic scintillator which is read 
out by wavelength-shifting (WLS) bars coupled to photomultipliers. Some of the features 
which make plastic scintillator attractive are: 

• High efficiency for crossing particles registration; 

• Fast scintillation timing; 

• Ease of calibration procedure; 

• Long-term stability and reliability; 

• Production potential - the plastic scintillator facility at IHEP, Protvino allows to 
produce the full amount of scintillator within 2-3 years; 

• Low maintenance — the plastic scintillator detector is quite robust and will require 
little maintenance in the underground tunnel conditions. 

The proposed detector consists of scintillating counters (10 mm thickness, 50 cm width 
and up to 6 m length) which cover the walls of the UNK tunnel. The proposed transverse 
coverage of the UNK tunnel is presented in Fig. 6. 

The adjacent scintillation counters are coupled with each other through the WLS- 
shifter (3 mm thickness, 10 mm width) for light collection along the 6 m length (see 
Fig. 7) and thus all counters in one 6 m cylindrical section form the continuous circular 
chain. 

Scintillation counters based on polystyrene with PTP and POPOP fluorescence dopants 
can be produced by extrusion technology at the IHEP, Protvino scintillator production 
facility. The PMMA-based WLS-shifter bars with Kumarin-30 Fluorescence doping will 
be used. The light attenuation length for such WLS-shifter is of 3 m (19], ^0|. The 
scintillator (decay time 2.3 ns; emission peak 420 nm; attenuation length 20 cm at 
420 nm) and WLS bars (decay time 2.7 ns; emission peak 460 nm; attenuation length 
up to 3 m at 460 nm) are wrapped in white material (e.g. TYVEK) and black paper. 
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Fig. 6. The proposed transverse coverage of the UNK tunnel. 




Shifter (thickness 3-5 mm) 
514 



n 
l 




I / Sc plate (L = 6m) j 



12° 



Fig. 7. The scintillation counters coupling. 

The light is collected from both sides of the WLS-shifter bar using the 1" PMTs. A green 
extended phototube FEU-115M from MELZ (Moscow, Russia) |21] can be used for light 
collection. The PMT has an average quantum efficiency of 15% at 500 nm with a gain 
around 10 6 . The estimations show that the signal for muon, crossing the scintillating strip 
near the PMT will be ~ 25 ph. el. Signals from phototubes are delivered to the QDCs 
and TDCs. 

The coordinate resolution is estimated as ~ 30 cm along the tile and ~ 10 cm in 
transversal direction using the TDC and QDC information. 

Each cylinder section with diameter 5.1 m and length 6 m, contains 32 counters. 
Trigger is defined as coincidence of at least two nonadjacent scintillation counters and 
performed only within neutrino source spill time. 

In the Table 4 main features of the UNK underground detector are given. 
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Table 4. Main parameters of the UNK detector. 



Item 


Value 


counter dimensions 


1U x oUO x oOOO mm 


time resolution 


1 ns 


coordinate resolution 


or „ „„ 

25 cm 


memciency 


4-10 


number of counters in the cylinder section 


on 
61 


total number of sections 


3500 


total number of counters 


112000 


total weight of scintillator 


3400 Ton 


total volume of the UNK tunnel 


400000 m 3 



4 Physics performance 
4.1 Statistics 

To estimate the statistics of the experiment we performed a Monte-Carlo simulation of the 
detector response for the simple oscillation approach of vacuum oscillations Piv^ — > v^) = 
1 — 7 M • sin 2 (1.27 • Am 2 ■ L/E). We used the coverage of the UNK tunnel shown in Fig. 6. 
We concentrated our study on the selection of and CC events. The following criteria 
were applied for selection: 

1. Presence of two hits in two non-adjacent scintillation counters. 

2. Time difference between hits greater than 10 ns. 

3. "Tracks" pointing to the neutrino source. 

4. Angle between neutrino and "track" directions less than 30° in each plane, and 
"track" direction below horizon. 

5. 0.7 < = v/c< 1. 

In the Table 5 the CC and NC event rates for different neutrino energy settings are 
presented. 

Table 5. Number of events in case of no oscillation for 10 7 s running time for each energy 
setting. 



E v , GeV 






CC 


NC 


CC 


NC 


GSI 


JHF 


GSI 


JHF 


GSI 


JHF 


GSI 


JHF 


1.7 


560 


160 


12 


4 


150 


50 


3 


0.7 


3.5 


3400 


1100 


70 


18 


1100 


360 


15 


3.3 


7.0 


12000 


4200 


200 


40 


2800 


1000 


30 


9 


14.0 


2800 


1000 


34 


9 


380 


130 


4 


1.1 
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As it is seen a reasonable number of CC events can be selected for the NBB with 
E v > 2 GeV within one year of running (10 7 s) for the JHF and GSI cases. 

Fig. 8 shows the dependence of the detector effective mass for the CC and NC events 
as a function of the neutrino energy. Fig. 9 shows (a) CC v^jv^ effective mass ratio and 
(b) z/ M NC/V M CC and z/^NC/z/^CC event ratios. 
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Fig. 8. Effective mass for (a) CC and events, (b) NC and events as a function 
of the neutrino energy. 




As it is seen the CC effective mass linearly increases with the increase of E v reflecting 
the muon range increase. Slight increase of the NC effective mass is seen as well. The 
z/^CC/z/^CC ratio for effective masses is ~ 1.3 due to the different y-dependence of z/^CC 
and z/^CC cross sections. The NC/CC ratio is about few percent for E u ~ 2 GeV. The 
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Table 6. Estimation of running time needed for 10 3 events at each energy setting. 



v source 


E v , GeV 






Nev 


pot,-10 2i 


time, TO 7 s 


Nev 


pot,-10 2i 


time,-10 7 s 


JHF 


5.7 


1000 


0.35 


0.35 


1000 


1.4 


1.4 


7.1 


1000 


0.23 


0.23 


1000 


1.1 


1.1 


9.4 


1000 


0.32 


0.32 


1000 


1.3 


1.3 


14.0 


1000 


1.00 


1.00 










Total 


1.9 


1.9 


Total 


3.8 


3.8 


GSI 


2.0 


1000 


0.35 


1.04 


400 


0.56 


1.67 


4.0 


1000 


0.07 


0.22 


1000 


0.24 


0.72 


8.0 


1000 


0.02 


0.07 


1000 


0.12 


0.36 




Total 


0.44 


1.33 


Total 


0.92 


2.75 



efficiency for z/ e CC and z/ T CC is comparable with that for z/^NC. Thus such kind of detector 
is a natural u^CC event selector. 

To estimate the experimental sensitivity of the oscillation pattern measurements we 
select three-four neutrino energy settings where oscillations are maximal or minimal for 
Am 2 = 2.5 • 10~ 3 eV 2 . The statistics for each setting was restricted to 1000 CC events 
assuming that there is no oscillation. This restriction is reasonable assuming systematic 
uncertainties of about ~ 3% for the expected number of events using the proposed NBB. In 
the Table 6 we present these energy settings with the expected statistics (no oscillation), 
values of pot's and exposure times for and cases for the JHF and the GSI. As 
it is seen from this table the oscillation pattern measurements can be carried out in a 
reasonable time even in the worse case. 

4.2 Background 

The experiment sensitivity depends on the background level. It is expected that the 
main background is due to the cosmic muons. Three sources of this background can be 
considered: 

• fibl cosmic muons which coincide in direction with accepted solid angle of muons from 

the neutrino source; 

• fib2 cosmic muons which interact in the surrounding soil and produce through-going 

secondary particles within the accepted neutrino source direction; 

• /i/33 cosmic muons which move in the direction opposite to the accepted solid angle for 

the neutrino source but with wrong TOF identification. 

The fxbl background can be estimated knowing the muon flux as a function of the zenith 
angle 9. Fig. 10 shows this muon flux for the UNK detector for an average depth of 50 m. 
It is seen that the flux sharply drops down to the value of 4 • 10~ 7 //•m -2 -sr -1 -s -1 at 
9 = 90°. This constant value for 9 > 89.5° is defined by muons from atmospheric neutrino 
interactions. It means that we should accept only muons below horizon (9 > 90°). This cut 
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20 40 60 80 100 120 140 160 180 

, degree 

Fig. 10. Muon flux as a function of zenith angle 9 for the UNK underground detector 
(50 m depth in average). 

reduces the statistics for the GSI neutrino source by ~ 20%. We estimate this background 
as N^ b i = T ■ ■ Q ■ S ■ D — 0.8 \i for an exposure time T = 10 7 s (one year), a solid 
angle Q = 1 sr, a muon flux = 4 ■ 10~ 7 /i-m _2 -sr _1 -s~ 1 , a detector area S = 10 5 m 2 
and a duty factor D = 2 • 10~ 6 . 

The [ih2 background is estimated using the MACRO background measurements |22 . 
It was found that each down-going muon produces the 2 • 10 -5 registered up-going par- 
ticles. It is estimated that for the UNK detector this background is N^2 = 40 /i for one 
year exposure (T = 10 7 s). Therefore, it is the most significant source of cosmic muon 
background. It can be further reduced using the RF structure of the beam. 

We estimate the fib3 background as the value of = 0.2 fi for 10 7 s assuming 

$(/x) = 1 /i-m~ 2 -sr~ 1 -s~ 1 and the Qa separation for the TOF system (10~ 7 probability for 
the Gaussian). 

For the first glance the cosmic muon background is not so severe for our experimental 
conditions and is < 20% for the worse case at low energies. However, it should be noted 
that the background is sensitive to the detector characteristics like tails of coordinate and 
time resolutions which were not taken into account. Obviously, a direct measurement of 
the background in the UNK tunnel is needed. 



4.3 Sensitivity 

To estimate the sensitivity of the UNK neutrino experiment to disappearance oscilla- 
tion parameters we used the values of Am 2 = 2.5 ■ 10~ 3 eV 2 and 1^ = 0.95 and selected 
the neutrino energy settings with the expected numbers of CC events (no oscillation) 
presented in Table 6. Fig. 11 shows the simulated probability values P = N osc /N exp with 
error bars, the fitted curve and the obtained ideal oscillation curve for the GSI and the 
JHF examples for the identification and measurement of the neutrino oscillation 

patterns. 
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Fig. 11. Simulated oscillation pattern (full points), fitted (solid line) and ideal (dashed 
line) oscillation curves for the GSI and the JHF cases using of data, presented in 
Table 6. 



The estimated sensitivity for the GSI case is 



and for the JHF case is 



<r Am 2 = 2.7 • 10 
<7j„ = 0.01, 



-5 



eV and 



O Am 2 



= 1.5-10' 
= 0.01. 



eV and 



So, in this experiment the relative errors ~ 1% can be achieved for Am 2 and 1^. 
These sensitivities can be reached in the Am 2 regions as presented in Table 7. 

Table 7. The Am 2 region accessible to the experiment. 



v source 


Am 2 in 10~ 3 eV 2 


GSI 


1.5 < Am 2 < 6 


JHF 


0.5 < Am 2 < 4 



The measurement of A M = P{y^ — > u^) — P{p^ — > v^) as a function of the neutrino 
energy can allow to estimate the values of effective CPT violation parameters 5^ = 
Am 2 — Am 2 and e M = 1^ — Ip with a relative precision of ~ 2% and of ~ 1%, respectively. 

Such a precision allows to observe oscillation or decay neutrino patterns, to measure 
precisely the neutrino oscillation parameters in the simple case of oscillation, to search for 
fake and genuine CPT violation effects at a few percent level and to search for complicated 
oscillation behaviour. 
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4.4 Terrestrial matter effects 



The analysis given before was based on a simple approximate formula Piv^ — > v^) = 
1 — Ifj, sin 2 (1.2 7 Am 2 • L/E) for two-neutrino oscillation in vacuum, which contain two in- 
dependent parameters only. In this section we consider more realistic approach including 
terrestrial matter effects. We restrict themselves only by analysis in the framework of 
mixing between three known neutrino. This scheme is compatible with the existing data 
on solar and atmospheric neutrino observation as well as reactor neutrino experiments. 
Four- neutrino scheme with light sterile neutrino is required only by LSND experiment, 
which is still not confirmed. Moreover, recent analysis tells us that this scheme is un- 
satisfactory as explanation of all experimental evidences for neutrino oscillations and is 
marginally acceptable [p3fl . 

In the case of three-neutrino scheme, unitary matrix which transforms neutrino mass 
eigenstates into flavour ones can be written in the standard form 

C12C13 S12C13 s* 3 

U = ( -S12C23 - C12S13S23 C12C23 - S12S13S23 C13S23 J ■ (1) 

S12S23 ~~ c 12§13 c 23 — C12S23 — S12S13C23 C13C23 

where §13 = s^e 16 , = sin^- and = cos 6^-. The matrix has four independent 
parameters: three mixing angles 9ij (i < j) and CP violating phase 5. In general there 
exist two additional Majorana specific phases, but they do not contribute to the lepton 
number conserving oscillations and are omitted in Eq. |I[ Four parameters in matrix 
U, together with two squared mass differences for neutrino, say Am^ = m\ — m\ and 
Am§ 2 = m 2 — m 2 , form six parameter set which defines full oscillation pattern. At present, 
the best fit values for these parameters lie in two distinct statistically significant regions 
called large mixing angle (LMA) and low mass (LOW) |23|]. They are the two solutions 
which we will keep in mind in further investigation. We do not intend to explore the all 
allowed parameter space. Rather, we are going to answer general questions such as the 
principal possibility to measure CP violation parameter 5 or the possibility to distinguish 
LMA and LOW solutions. 

Differential equations which describe propagation of neutrino flavour states on a dis- 
tance L look as follows 



._d_ 
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-Am 2 ! \ / A(L) 

0^ C/t + 
Am 2 2 y \ 









)] 


1 











, (2) 



A(L) = -L Gf N aP (L). 

Here E is the energy of neutrino beam, A(L) is the potential induced by coherent in- 
teraction of electron neutrinos with electrons in the matter, Gf is the Fermi constant, 
Na is Avogadro's number and p(L) denotes the Earth matter density. Equations for 
antineutrinos are obtained from Eq. ^| by replacement (U,A) — > (U*,—A). Fake CPT 
violation appears due to the different sign at A for neutrino and antineutrino. In the 
vacuum P(z/ M -> v^) = P(p^ -> v^. 

In our analysis we solve these equations numerically with density profile p(L) taken 
from the Preliminary Reference Earth Model ||24j| . As for the mixing angles and neutrino 
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squared mass differences, we choose two sets of the parameters corresponding to the 
aforementioned LMA and LOW regions 



Am 2 21 = 1.4 ■ 1(T 4 eV 2 



|Am 2 2 | = 2.5 ■ 1(T 3 eV 2 



12 = 35° Point 1 (LMA) 



23 = 40° (3) 
013 = 13° 

Am§! = 1.0 • 1(T 7 eV 2 Point 2 (LOW) 
9 12 = 35°. 

These values will be used as reference points in parameter space but we will also consider 
variations in some of parameters in the experimentally allowed limits. 

To characterize the effects of fake CPT violation due to the non-zero matter density 
it is convenient to introduce asymmetry between - oscillation probabilities 

A = P K -> V A ~ Zl V A (a) 
CPT p(^-„ M )+p(^-p M )- [) 

Before discussion the results of our calculations we make some general remarks. 

1. Without loss of generality one may take for mixing angles and phase in matrix U the 
following values 

7T 

0^ G [0, — ] and 5 G [ — 7r, 7r] . 

From Eq. |2| one can show that, in the constant matter density approximation, P(z/ At (z/ /i ) — > 
Vnipfj)) does not alter under the change 5 — > —5. The experiment under discussion could 
not distinguish 5 and —5. 

2. When any of the mixing angles in matrix U is zero, the dependence on phase 5 in 
oscillation probabilities disappear. In the experimentally allowed parameter space is 
the smallest angle and it is the only one which is compatible with zero (0 < s 2 3 < 0.05 
and zero value minimizes the x 2 function for global fit of present data 

3. In the limit Am|i — > the dependence on angle 9±2 and phase 5 is dropped out in Eq. ||. 
Therefore, there is no chance to determine these two parameters in discussed experiments 
when LOW solution is realized in Nature. 

4. In the limit Am^ — > 0, the changing in sign of Am| 2 is equivalent to passing from 
Eq. ^ for neutrinos to equations for antineutrinos, i.e., P{v^ — > Am^) = Pip^ — > 

-Amy. Therefore, the sign of asymmetry Aqpt, if not zero, may reflect the sign of 
Amg 2 even in the case Am 21 7^ 0. 

5. For the LOW solution (for which Am 21 = is very good approximation) the matter 
influence on oscillation pattern is sensitive to the value of i3 . At 9 13 = Eq. || decouples 
into two pieces: electron neutrino does not take part in oscillations while muon and tau 
neutrinos oscillate as in vacuum and Acpt vanishes. This is not so for the LMA solution, 
zero value of #13 does not lead to the decoupling due to the non-negligible value of Am 21 
and deviation of Acpt from zero may be significant. Hence, the measurement of near 
zero value of Acpt at all energies will point out on the LOW solution with small 6*13. 
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Numerical results In all figures presented below we use the oscillation probability 
smeared with energy, assuming NB like beam with energy spread SE/E = 0.15. 

Figs. 12 and 13 illustrate the matter influence on the oscillation profile for the case of 
GSI and JHF, respectively. All curves correspond to phase 5 = 0° and Am| > 0. The 
difference between neutrino and antineutrino oscillation curves is seen directly, especially 
it is impressive for the case of JHF. It is interesting to note, that in the last case there is 
not only difference in amplitude but also a phase shift, which reaches up to ~ 1 GeV. 

Asymmetry Aqpt is plotted in the Figs. 14 and 15 as function of beam peak energy. 
One can see that asymmetry has absolute maximum both for GSI and JHF cases. As for 
the JHF case, maximum in asymmetry appears near the position of the first maximum of 
the oscillation curve where the difference between neutrino and antineutrino probabilities 
is enormous. In the GSI case the absolute differences between neutrino and antineutrino 
at first minimum and maximum of oscillation curve are approximately the same and, 
therefore, maximum in asymmetry occurs at first oscillation minimum. 

Thick lines in Figs. 14 and 15 corresponds to positive Am| 2 , while thin lines to negative 
one. Solid lines in Figs. 14(b) and 15(b), which are for values of parameters at Point 2 
(Eq. [|), demonstrate the reflection of asymmetry under the changing in sign of Am^. 
Such a behaviour is kept, in the energy region near the asymmetry maximum, even for 
the parameters at Point 1, in spite of the not so large hierarchy between Am^ and Am^. 
Clear distinction between the positive and negative signs of Am| 2 is lost only at lower 
energies. One may conclude that at energies near the asymmetry maximum the sign of 
Am\ 2 is defined unambiguously, whatever solution, LMA or LOW, is realized in Nature. 

To illustrate the sensitivity to the angle #13 for LOW region we present in Figs. 14(b) 
and 15(b) the additional curves for asymmetry corresponded to sf 3 = 0.01 and 0.03. It is 
clear that the JHF, in contrast to the GSI case, is in a position to resolve small values of 
angle 9 13 . 

As was mentioned early, the experiment may be sensitive to the CP violating phase 
5 only for LMA solution. In Figs. 14(a) and 15(a) solid lines present the asymmetry 
for 5 = 0, while dashed lines are for 5 = 180°. One can see that at low energies, say, 
at ~ 1 GeV, the values of asymmetry are large enough and clearly separated for these 
boundary values of 5. It is in contrast with the behaviour of asymmetry for LOW region of 
parameters, where it approaches to zero at low energies. As a consequence, the observation 
of non-zero value of asymmetry at low energy will discriminate between LMA and LOW 
solutions, on the one hand, and will allow to determine the value of 5, on the other hand. 
The measurement of near zero value of Aqpt at low energies will be less informative: there 
exist some values of parameters when oscillation patterns for LOW and LMA solutions are 
very close. Fig. 16 illustrates possible worse case to study for the parameters of Point 1 
with 5 = 90° and Point 2. 

One should note that the using of NB like beam is important, first of all, for the precise 
measurement of the shape of oscillation curve and, hence, for the precise measurement of 
oscillation parameters like Am 2 or amplitude 7 M . In contrast with this, the low energy 
measurement does not require NBB. Rather, one should use WB like beam to compensate 
the decreasing in statistics at low energies. 
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LMA, GSI LOW, GSI 




1. 5. .1e2 1. 5. .1e2 

E, GeV E, GeV 

(a) (b) 

Fig. 12. Survival probability as function of beam energy for the GSI (L = 2000 km): 
(a) Point 1, (b) Point 2 (Eq. [3]). Thick solid line shows Piy^ — > z/ /t ) in matter smeared with 
8E V /E V = 0.15, thick dashed line shows smeared -P(^ — > in matter. For comparison, 
smeared vacuum oscillation curve (thin solid line) and vacuum oscillation curve without 
smearing (thin dashed line) are also plotted. All curves correspond to 5 = 0°. 



LMA, JHF LOW, JHF 




Fig. 13. The same as in Fig. 12, but for the JHF case (L = 7000 km). 
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LMA, GSI 



LOW, GSI 



0.2 z 




Fig. 14. Neutrino-antineutrino asymmetry as function of beam peak energy for the GSI 
(L = 2000 km): (a) Point 1, (b) Point 2 (Eq. |3|). Thick lines correspond to positive 
Am| 2 while thin lines to negative ones. Solid lines are for 5 = 0°, dashed lines are for 
5 = 180°. Wide grey strip shows errors in asymmetry with 5 = 0° and Am| 2 > 0. Errors 
correspond to statistics for 10 7 s running time. Dotted lines in (b) are for sf 3 = 0.01 and 
0.03 (Am^ 2 > 0); correspondence is as follows: the smaller the angle #13 the smaller the 
asymmetry. 



LMA, JHF LOW, JHF 




J7~ 5. .1e2 1. 5. .1e2 

E, GeV E, GeV 



(a) (b) 
Fig. 15. The same as in Fig. 14, but for the JHF case (L = 7000 km). 
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5. .1e2 1. 5. .1e2 

E, GeV E, GeV 



(a) (b) 

Fig. 16. Worse case to study: Point 1 and 5 = 90° vs. Point 2 and any 5. (a) -P(V M — > Vp)\ 
(b) Asymmetry. Two thick lines correspond to GSI [L = 2000 km), Points 1 and 2, two 
thin lines correspond to JHF (L = 7000 km), Points 1 and 2. 



LMA, GSI 



LMA, JHF 
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17. 

,2 



P(u„ 



v., 



vs. 



Pivp -> z^) at E 



and s^ 3 for Point 1: 



Am 2 21 



2.7-10" 



to s 



/; , . , ■ ;; ■ ,, , .... ^ 1 GeV in dependence on phase (5 
» for GSI, (b) for JHF. Regions for Am? 21 = 1.4 • 10" 4 eV 2 and 
4 eV 2 are marked by "1.4" and "2.7", respectively. One point corresponds 
(black circle) and five points with 5 = 0°, 45°, 90°, 135° and 180° (marked by 
the same symbols) are plotted for the values of sf 3 = 0.01 (boxes), 0.03 (circles) and 0.05 
(diamonds). Straight line grid presents different values of asymmetry Aqpt = (solid 
line), ±0.04, ±0.08, . . . (dashed lines). The errors for point sf 3 = 0.03 and 5 = 45°, which 
correspond to the statistics of 10 3 events in the absence of oscillations, are plotted. 
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The sensitivity to 5 depends on the value of s 2 3 and the maximal sensitivity is for 
maximal allowed value of sf 3 = 0.05. To comprehend the influence of sf 3 we calculated 
P{Vfx — * Vfj) and Pip^ — > V^) at beam energy E = 1 GeV for various values of s 2 3 and S. 
Fig. 17 presents the results of these calculations. It contains two distinct regions which 
correspond to Am^ = 1.4- 10~ 4 eV 2 and Am^ = 2.7- 10~ 4 eV 2 . One can see that for the 
case of GSI there is one-to-one correspondence between points in the plane (P(i/ /1 ), Pip^) 
and values (s 2 3 , 8). The typical errors, plotted for one point (s 2 3 = 0.03, 5 = 45°), allows to 
conclude that parameters s 2 3 and 5 could be, in principle, disentangled and determined 
simultaneously. The situation for the case of JHF is more involved. The one-to-one 
correspondence is lost and the relative variations in oscillation probabilities with s 2 3 and 
5 are smaller than in the GSI case. Qualitatively this can be understood as the increasing 
role of matter effects. The last ones become dominant even at relatively low energies, 
where the sensitivity to 5 is maximal. 

Considered in this section topics allow to make a comparison between GSI and JHF 
as the sites for neutrino beam source. If LMA solution is realized in Nature, GSI can 
provide more information on oscillation parameters in the low energy measurements. JHF 
is more sensitive to the parameters of LOW region at higher energies where the matter 
effects maximize asymmetry Aqpt- 

5 Conclusion 

The VLBL neutrino oscillation experiment (baseline of 2000-7000 km) with the UNK 
~ 1 MTon underground detector will provide: 

• Observation of the oscillation patterns for the and disappearance in the 0.5 • 
10~ 3 < Am 2 < 6 ■ 10~ 3 eV 2 range using the NB neutrino beam. 

• Direct observation of the matter effects by comparison of and oscillation curves. 
It allows to extract important physical information such as 

1. Measurement of the Am§ 2 and the intensity of oscillations with the accuracy at the 
level of ~ 1%. 

2. Measurement of sin 2 #i 3 down to the value of ~ 0.01. 

3. Unambiguous determination of the sign of Am 32 by the sign of asymmetry A C pt 
for sin 2 13 > 0.01. 

4. Distinction between LMA and LOW solutions outside the case of some exceptional 
values of mixing matrix parameters. 

5. Determination of CP violating phase 5 in the case of LMA solution. 

Thus, the experiment allows to prove the standard scenario of neutrino oscillations, to 
measure precisely essential oscillation parameters as well as to search for new phenomena. 
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